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Reliable predictions of light charged particle production in spallation reactions are important to correctly assess gas production in 
spallation targets. In particular, the helium production yield is important for assessing damage in the window separating the accel- 
erator vacuum from a spallation target, and tritium is a major contributor to the target radioactivity. Up to now, the models available 
in the MCNPX transport code, including the widely used default option Bertini-Dresner and the INCL4.2-ABLA combination of 
models, were not able to correctly predict light charged particle yields. The work done recently on both the intranuclear cascade 
model INCL4, in which cluster emission through a coalescence process has been introduced, and on the de-excitation model ABLA 
allows correcting these deficiencies. This paper shows that the coalescence emission plays an important role in the tritium and 3 He 
production and that the combination of the newly developed versions of the codes, INCL4.5-ABLA07, now lead to good predic- 
tions of both helium and tritium cross sections over a wide incident energy range. Comparisons with other available models are 
also presented. 
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1. Introduction 

Spallation reactions induced by high-intensity proton beams 
in a heavy-metal target are used to produce intense neu- 
tron fluxes for various applications such as spallation neutron 
sources, accelerator-driven sub-critical reactors, or radioactive 
ion beam facilities. These reactions produce large quantities of 
light charged particles (hydrogen and helium isotopes) which 
are a concern in spallation target design. Indeed, build up of 
gases (in particular helium) can lead to swelling and embrit- 
tlement of the window separating the accelerator vacuum and 
the spallation target and of other structural materials; tritium, 
a twelve-year half-life beta emitter, is a concern for radiopro- 
tection, especially in the case of liquid targets from which it 
can escape easily. It is therefore important that nuclear physics 
models implemented into the high-energy transport codes that 
are used by spallation target designers be able to reliably predict 
light charged particle production yields. 

Up to now, light composite particle production was poorly 
predicted by nuclear models (generally a combination of an in- 
tranuclear cascade model followed by a de-excitation model) 
implemented into high-energy transport codes QjJ, |2|, [3j |4J] . In 
Ref. jlj] for instance, the cross sections for tritium production 
on iron and lead, as a function of the incident proton energy, 
were compared with MCNPX o] calculations using three dif- 
ferent models: Bertini-Dresner Q] , which is the default op- 
tion, CEM2k from Ref. H and Isabel H] -Dresner. It was found 
that these models were not able to account for tritium produc- 
tion, except Isabel-Dresner but only for lead. Discrepancies 
up to a factor of five were observed. Also, the dependence of 



the cross section with the incident energy is often not properly 
reproduced. In Ref. J2l, experimental data from |Kj IT] on 
helium and hydrogen production cross sections as a function 
of the target charge, at different energies, were compared with 
FLUKA 11211 . A systematic underprediction by the model, more 
important for helium and heavy targets, was reported. In J3l, 
a systematic comparison of many different models with alpha 
production cross section from tantalum, tungsten and gold was 
performed, which also pointed out large discrepancies between 
the calculation results and the experimental data. 

During the last ten years, in Europe, an important effort has 
been devoted to the collection of high quality experimental data 
and, simultaneously, to the development of improved spallation 
models lfl3ll . As a result of this work, the combination of the 
intranuclear cascade model, INCL4.2 [14], and de-excitation 
model, ABLA [15], has been developed, tested against a large 
set of experimental data and shown to give globally better pre- 
dictions than models used by default (Bertini-Dresner |@, Q]) 
in high-energy transport codes such as MCNPX H. INCL4.2- 
ABLA is now available in MCNPX and GEANT4 With 
these models, the situation regarding the predictions of neutron 
emission, heavy evaporation and fission residues production 
could be considered as having been largely improved. However, 
important deficiencies were still remaining when light charged 
particles are concerned. For helium out of iron, a systematic un- 
derprediction by a factor 2 to 3 by INCL4-ABLA was reported 
in Ref. |3[]. Since in ABLA only nucleons and alphas can be 
evaporated, INCL4-ABLA does not produce any tritium. This 
is obviously a major shortcoming since, for instance in calcula- 
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tions of radioactivity of the MEGAPIE target 11711 . tritium has 
been identified as a major contributor around 10 years after ir- 
radiation when using models as Bertini-Dresner or CEM. 

New versionsof INCL4 and ABLA (INCL4.5 0, and 
ABLA07 iT " 



2 1H respectively) have recently been released, 
which, among other improvements, were developed to remedy 
this situation: cluster emission, already introduced in the in- 
tranuclear cascade model in Ref [22] to account for the high- 
energy tail observed in the experiments, has been refined; evap- 
oration of deuterium, tritium and 3 He is now taken into account 
in ABLA together with the possibility of emitting intermediate 
mass fragments, and barriers for light charged particle emission 
have been established on better physics basis. In this paper, we 
present results concerning improved predictions of helium and 
tritium obtained with these new versions. The work on the mod- 
els has been done in the framework of the NUDATRA domain 
of the EUROTRANS European FP6 project H, whose objec- 
tive is to provide improved simulation tools for the design of 
ADS transmuters. 

2. Cluster emission in INCL4 

In light composite particle energy spectra, two different com- 
ponents are generally observed, as shown by the experimental 
data from Herbach et al. Ill ill plotted in Fig. Q~| a low-energy 
isotropic peak coming from the evaporation stage, and a high- 
energy forward-peaked tail. An attempt to separate the two 
components was done by the authors of Ref. Ill ill . They ex- 
tracted the relative contribution of this high-energy component, 
called pre-equilibrium in their paper, to the total yields of the 
different light composite particles, as a function of the target 
charge, at 1.2 GeV. This showed that, while 4 He is produced 
predominantly by evaporation, for the other light composite 
particles the so-called pre-equilibrium contribution is far from 
being negligible, reaching even 60% for 3 He on high-Z targets. 

The level of this contribution in the production cross sec- 
tions of tritium and 3 He shows how important it is to be able 
to account for the high-energy tail with the models. This is the 
reason why a coalescence mechanism leading to light compos- 
ite particle emission has been implemented into INCL4, first in 
Boudard et al. B22ll . and then modified in order toget a better ra- 
tio between the different species of clusters jl8lll9ll . The cluster 
emission mechanism is based on surface coalescence in phase 
space, i.e. on the assumption that a cascade nucleon ready to 
escape at the nuclear surface can coalesce with other nucleons 
close enough in phase space and form a cluster. The param- 
eters of the model include the volume of the phase space cell 
in which nucleons should be to form a cluster and the distance 
from the surface at which the clusters are built. All possible 
clusters up to a given mass number are formed and the priority 
is given to the one with the lowest excitation energy per nu- 
cleon. The selected cluster is emitted only if it succeeds to go 
through the Coulomb barrier. In INCL4.5, used in this paper, 
clusters up to mass 8 are considered. 

The other modifications compared to the original INCL4 
model M22ll involves the introduction of a potential for pions 
and of an isospin and energy dependent nuclear potential iTlsIi . 
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Figure 1: Light composite particle double-differential cross sections, in the 
p+Ta reaction at 1.2 GeV, measured by Hill and compared with the results 
of INCL4.5-ABLA07. The cross sections are mutliplied by 10°, KM and 10 2 , 
successively, starting from the smallest angle. The contribution due to the coa- 
lescence process in INCL4.5 is given by the dashed line. 



Despite the fact that intranuclear cascade models are based on 
assumptions that are in principle not valid below 100-150 MeV, 
INC models are nevertheless used in transport codes below 
these energies, for instance when evaluated data files for the 
considered nuclei do not exist. Therefore, it is important that 
the models give reasonable results also at low energies. In the 
case of INCL4.5, in order to get better results at low energies, 
in particular concerning the total reaction cross section, special 
attention has been paid to the treatment of the first nucleon- 
nucleon collision and Coulomb distortion in the entrance chan- 
nel has been introduced 10. 

In Fig. Q] results on light composite particle double- 
differential cross sections obtained with INCL4.5 coupled to 
the de-excitation model ABLA07 ll20ll are compared with the 
experimental data from [11] measured at 1.2 GeV in the p+Ta 
reaction. The dashed curve shows the contribution coming from 
the coalescence process in INCL4.5. It can be seen that the 
high energy tail is very well reproduced by this process for 
all species of light charged particles, although a slight over- 
prediction of 3 He and 4 He at the highest particle energies can 
be observed. The angular dependence is also reasonably repro- 
duced. It can be noticed, as already realized by the authors of 
Ref. JU, 

that the cascade contribution to the total cross section 
(dashed curve) is very important in the case of tritium and 3 He 
contrary to the case of 4 He. A similar agreement, not shown 
here, is obtained when comparing the model with p+Au data 
from Ref. ll24ll at 2.5 GeV (this and further results can be found 
in the benchmark of spallation models recently organized by 
IAEA Jll). 



2 




20 40 60 

E (MeV) 



20 40 60 

E (MeV) 



Figure 2: Same as Fig.fT]but for the p+Fe reaction at 62 MeV, measured by 1 26]. 

Fig. [2] gives an example of the results obtained at low inci- 
dent energies, here 62 MeV for the system p+Fe compared to 
the data from j26ll . The agreement is not as good as at high in- 
cident energies, in particular at forward angles where the slopes 
of the calculated energy spectra are too steep, but still satisfac- 
tory. At this incident energy, the cascade coalescence contribu- 
tion is clearly dominant except at backward angles for d, t, 3 He 
and for A He at all angles. It should be stressed that the good 
agreement achieved for different target masses on a wide range 
of incident energies is obtained with the same set of parameters 
in the coalescence model, chosen once for all. 

3. Recent improvements of ABLA 

Simultaneously to the work on INCL4, the authors of 
ABLA have produced a new improved version of the model, 
ABLA07 1 20]. This version now allows evaporation of all the 
types of light charged particles from p to 4 He but also of in- 
termediate mass fragments. It uses improved parameterizations 
of inverse reaction cross sections and Coulomb barriers in or- 
der to better reproduce experimental particle energy spectra. 
A simultaneous break-up (multifragmentation) mechanism has 
been added for systems overcoming a certain limiting tempera- 
ture 12711 . The fission part has also been modified. More details 
as well as examples of the improved agreement with experi- 
mental data can be found in Ref. |2~lll . 

It can be seen in Fig. Q] that the low energy part of the deu- 
terium and tritium spectra, corresponding to the evaporation 
peak, is very well reproduced by ABLA07. For 3 He, the situa- 
tion is less good since the production yield is a little too small 
and the evaporation peak is not broad enough, especially on the 
low-energy side. 4 He is intermediate between tritium and 3 He. 



Since INCL4.5-ABLA07 also agrees very well with the neu- 
tron evaporation spectra measured on Pb at 1.2 GeV (results 
not shown here but which can be found in the IAEA bench- 



mark [25]), we are rather confident that excitation energy distri- 
bution at the end of the cascade stage is correct. Therefore, the 
slight discrepancy observed for helium isotopes is more likely 
an indication that the barriers and/or tunneling through the bar- 
rier in case of the helium emission in ABLA07 still have to be 
improved. 

4. Total tritium production yields 

Figs. [3] and [4] present the results obtained with INCL4.5- 
ABLA07 (solid red line) for the total production yields of tri- 
tium on iron and lead targets, respectively, as a function of the 
incident proton energy, compared to the available experimen- 
tal data SIIi[3S|3l[3l[3l[3l,[iIl[il- Calculations with 
other models implemented in MCNPX, namely Bertini-Dresner 
(blue line), which is the default option, and the new version of 
the CEM model (green line) from S. Mashnik, CEM03 Q- 
For both iron and lead, INCL4.5-ABLA07 gives a good agree- 
ment with the data all along the energy range. This combination 
of models is definitely better than Bertini-Dresner which has a 
totally wrong energy dependance and even than CEM03 which 
largely overestimates tritium production at high incident ener- 
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Figure 3: Tritium production cross sections in iron calculated with the new ver- 
sions INCL4.5-ABLA07 (solid red line) compared to data measured by differ- 
ent groups (Refs. [28, 29, 30, 3ll[ll[ll[Tl) and to calculations with Bertini- 
Dresner (blue) and CEM03 (green) models. The dashed, resp. dashed-dotted, 
lines show the respective contributions from the cascade coalescence and evap- 
oration processes. 

The respective contributions of the cascade coalescence 
(dashed line) and evaporation processes (dashed-dotted line) in 
the tritium production cross sections are also shown in the fig- 
ures. It can be seen that, at low energies, tritium is produced 
nearly exclusively by the coalescence process, evaporation be- 
coming dominant only above 600 MeV for Fe and 1.7 GeV 
for Pb. This illustrates the importance of having models able 
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Figure 4: Same as Fig.[3]but for lead. Data from Refs. |28, 29, 3fJ[ll[Tl[3l 



to account for the non-evaporative component of the spectra 
and explains why models such as Bertini-Dresner, or ISABEL- 
Dresner not discussed in this paper, cannot give good predic- 
tions, especially at low incident energies. In CEM03, high- 
energy clusters are produced either in a coalescence process or 
in the pre-equilibrium stage. 

5. Helium production yields 

Fig.|5]shows the result for the 3 He and 4 He production cross 
section in iron as a function of the proton incident energy calcu- 
lated with INCL4.5-ABLA07, with the respective contributions 
of cascade coalescence and evaporation, compared to the exper- 
imental data from [iil|35|]. Like tritium, 3 He is predominantly 
originating from the cascade coalescence process up to about 1 
GeV. On the other hand, 4 He is produced essentially by evap- 
oration, the cascade contribution representing at most 20% at 
the lowest energies. This illustrates the difference in the two 
physical mechanisms: in the evaporation process, the emission 
of tightly bound particles (alphas) is favored compared to more 
loosely bound ones (d, t, 3 He), while the coalescence is mainly 
governed by the number of nucleons in a cluster (there are more 
clusters of mass 2 than 3, 3 than 4), the binding energy of the 
clusters playing only a minor role. 

The combination of the two models reproduce very well the 
3 He cross section on the full energy range while it underpre- 
dicts 4 He a little bit. Since 4 He is predominantly produced by 
evaporation, this underestimation could be related with the fact, 
mentioned in section[3] that ABLA07 gives a too narrow evap- 
oration peak for helium isotopes. However, this observation 
was made for Ta and the situation could be different for a much 
lighter nucleus like Fe. Unfortunately, the double-differential 
cross sections on Fe measured by the NESSI group are not 
available. In Fig. [6] the results given by the CEM03 code are 
also compared to the experimental data and to our model. It 
can be seen that, at high incident energies, CEM03 overpre- 



Figure 5: 3 He (blue) and 4 He (red) production cross sections in iron calcu- 
lated with INCL4.5-ABLA07 compared to experimental data from (Til Si- 
The dashed, resp. dashed-dotted, curves give the contribution from the cascade 
coalescence, resp. evaporation emission. 



diets 3 He, as it was the case for tritium, but is slightly better 
than INCL4.5-ABLA07 for 4 He. 
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Figure 6: 3 He (dashed) and 4 He (solid) production cross sections in iron cal- 
culated with INCL4.5-ABLA07 (red) compared to data measured by different 
groups (Refs. |Tj,[33]) and t0 calculations with CEM03 ( green). 



Since what matters for material damage assessment is the to- 
tal helium production yield, the helium production cross section 
is shown in Fig. [7] for different models compared to the avail- 
able experimental data. The first remark is that the new version, 
INCL4.5-ABLA07, represents a clear improvement compared 
to the previous one, INCL4.2-ABLA, and to Bertini-Dresner, 
which, being the default option of MCNPX, is used by most 
of the users of the code. In INCL4.2-ABLA the step observed 
around 100 MeV was due to a forced absorption process added 
at low incident energies in order to obtain the right total reac- 
tion cross section. In the new model this is no longer necessary 
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since it naturally gives the correct total reaction cross section, 
as mentioned in section [2] As already observed, CEM03 pre- 
dicts somewhat larger cross sections above 500 MeV than our 
model. In view of the discrepancy between the different sets of 
experimental data, it is difficult to conclude which model is the 
best. 
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Figure 7: Helium production cross sections in iron calculated with INCL4.5- 
ABLA07 (red) compared to data measured by I'll, 35] and to calculations with 
INCL4.2-ABLA (dashed pink), CEM03 (green) and Bertini-Dresner (blue) 
models. 
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Figure 8: 4 He production cross sections in lead calculated with INCL4.5- 
ABLA07 (red) compared to data measured by 1371 11 lL 1381 13911 and to calcu- 
lations with INCL4-ABLA (dashed pink), CEM03 (green) and Bertini-Dresner 
(blue) models. 

Regarding 4 He production on lead, Fig. [8] shows the com 
)arison of our model with the experimental data from J37 , 11 
As for iron, INCL4.5-ABLA07 and 



38, 39] and other models. 
CEM03 give the best agreement with the experimental data, 
CEM03 producing a little more A He than our model at high 
incident energies. Bertini-Dresner has a much too steep depen- 



dence with incident energy, leading to a severe underestimation 
of helium at low energies. Compared to the INCL4.2-ABLA 
version, our new model gives similar results except at low in- 
cident energies where the step around 100 MeV has been sup- 
pressed. 
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Figure 9: Cumulative 3 He production cross sections in lead calculated with 
INCL4.5-ABLA07 (red) compared to data measured by f37|.[Tl|.[38l[39|l and to 
calculations with CEM03 (green). 

There do not exist many data on the production of direct 3 He 
from lead. Therefore, we have compared the models to the mea- 
surement by Leya et al. [38] of cumulative yields of 3 He, i.e. 
direct 3 He plus 3 He from the decay of tritium, as a function of 
incident energies. The results are shown in Fig.|9]in which the 
calculated cross sections by INCL4.5-ABLA07 and CEM03 of 
3 He and tritium have been added up to be compared to the ex- 
perimental data. The same was done for the Herbach data point 
displayed in the figure. The same conclusions as for tritium 
production in Fig. [4] can be drawn, i.e. that our model gives 
a better agreement to the data than CEM03. This actually is 
not surprising since tritium represents around 80% of the cu- 
mulative yield. This observation is rather an indication that the 
different sets of experimental data are consistent. 

6. Conclusion 

In this paper, we have compared the predictions of the 
newly developed versions of the INCL4 and ABLA mod- 
els, INCL4.5-ABLA07, to the helium and tritium produc- 
tion double-differential cross sections and excitation functions 
found in the literature. It has been shown that the combination 
of the new versions represent a definite improvement compared 
to the former models, which did not evaporate tritium and 3 He 
and were unable to account for the high energy tail of the par- 
ticle spectra. In particular, it was found that the coalescence 
process, added in the cascade model to describe the emission of 
high energy clusters, is very important to get correct production 
cross sections of tritium and 3 He since it represents the major 
part of the cross section up to rather high energies. 
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Comparisons with the other models available in MCNPX, 
Bertini-Dresner and CEM03, have also been shown, which def- 
initely rules out the first one although it is the most often, 
blindly, used (default) option. CEM03 was found to be less 
good than our model on tritium and 3 He production, especially 
at incident energies above 300-400 MeV. For A He or total he- 
lium cross sections, it is difficult to conclude which of the two 
models is the best in view of the discrepancies between the dif- 
ferent sets of data. 

The INCL4.5-ABLA07 combination of models will be soon 
available in MCNPX [40]. In view of its validation over the 
wide energy range presented in this paper, this should allow 
reliable simulations of, for instance, tritium build up in liquid 
metal spallation targets or gas production in material structures. 
Also, the good prediction of the high-energy tail of the emitted 
light-charged particles, should guaranty a reasonable assess- 
ment of the secondary reactions potentially induced by these 
particles. 

Acknowledgments 

This work was partly supported by the FP6 Euratom project 
EUROTRANS/NUDATRA, EC contract number FI6W-CT- 
2004-516520. 



References 



[i] 



[2] 

[3] 
[4] 

[51 
[61 
[71 
[81 
[9] 
[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 



[18] 



[19] 



B. Rapp, J.C. David , V. Blideanu, D. Dor, D. Ridikas, N. Thiolliere. 
Proceedings of International Workshop on Shielding Aspects of Acceler- 
ators, Targets and Irradiation Facilities (SATIF-8), Pohang, South Korea. 
May 22- 24, 2006 (2006); EURISOL DS/Task5/TN-06-04 (2006). 
L. Pienkowski, W. Gawlikowicz, D. Hilscher, EURISOL internal report, 
EURISOL DS/Task5/TN-06-05 (2006). 

S. Leray et al, Nucl. Instrum. Methods A 562 (2) (2006) 806. 

C. H.M. Broders and A.Yu. Konobeev, Nucl. Instrum. Methods B 234 
(2005) 387. 

J. S. Hendricks et al, AIP Conf. Proceedings 768 (2005) 1188. 
H. W. Bertini, Phys. Rev. 131 (1963) 1801. 

L. Dresner, Oak Ridge report. Technical report, ORNL-TM-196, 1962. 
S.G. Mashnik et at, Advances in Space Research 34 (2004) 1288. 
Y. Yariv and Z. Fraenkel, Phys. Rev. C 20 (1979) 2227. 

D. Hilscher et al, Journal of Nuclear Materials 296 (2001) 83. 
CM. Herbach et al, Nucl. Phys. A 765 (2006) 426. 

Fasso, A., Ferrari, A., Ranft, J. and Sala, PR., 2005, FLUKA: a multi- 
particle transport code, CERN-2005-10, INFN/TC.05/11. 
HINDAS Final Report, EU Contract FIKW-CT-0003 1 , J.P. Meulders, A. 
Koning and S. Leray ed. (2005). 

A. Boudard, J. Cugnon, S. Leray, C. Volant, Phys. Rev. C 66 (2002) 
044615. 

A.R. Junghans, M. de Jong, H.-G. Clerc, A.V. Ignatyuk, G.A. Kudyaev, 
K.-H. Schmidt, Nucl. Phys. A 629 (1998) 635. 

S. Agostini et al., GEANT 4-a simulation toolkit, Nucl. Instr. and Meth. 
A 506 (2003) 250. 

S. Lemaire et al., Proceedings of the International Conference on Nuclear 
Data for Science and Technology, April 22-27, 2007, Nice, France, Ed. O. 
Bersillon, F. Gunsing, E. Bauge, R. Jacqmin, and S. Leray, EDP Sciences, 
2008, p304. 

A. Boudard and J. Cugnon, Proceedings of the Joint ICTP- 
IAEA Advanced Workshop on Model Codes for Spallation Re- 
actions, ICTP Trieste, Italy, 4-8 February 2008, D. Filges, S. 
Leray, Y. Yariv, A. Mengoni, A. Stanculescu, and G. Mank 
Eds., IAEA IN DC(NDS)-5 30, Vienna, August 2008, p29. 
http://www-nds.iaea.org/reports-new/indc-reports/indc-nds/indc-nds-0530.pdf 
J. Cugnon, A. Boudard, S. Leray, D. Mancusi, Proceedings of the Interna- 
tional Topical Meeting on Nuclear Research Applications and Utilization 
of Accelerators, 4-8 May 2009, Vienna. 



[20] 



[21] 



[22] 

[23] 

[24] 
[25] 
[26] 
[27] 

[28] 

[29] 
[30] 
[31] 
[32] 



[33] 
[34] 
[35] 

[36] 
[37] 
[38] 

[39] 
[40] 



A. Kelic, M.V Ricciardi and K.-H. Schmidt, Proceedings of the 
Joint ICTP-IAEA Advanced Workshop on Model Codes for Spal- 
lation Reactions, ICTP Trieste, Italy, 4-8 February 2008, D. 
Filges, S. Leray, Y. Yariv, A. Mengoni, A. Stanculescu, and G. 
Mank Eds., IAEA INDC(NDS)-53 0, Vienna, August 2008, p!81. 
http://www-nds.iaea.org/reports-new/indc-reports/indc-nds/indc-nds-0530.pdf 
M.V. Ricciardi, A. Kelic, K.-H. Schmidt, Proceedings of the International 
Topical Meeting on Nuclear Research Applications and Utilization of Ac- 
celerators, 4-8 May 2009, Vienna. 

A. Boudard, J. Cugnon, S. Leray, C. Volant, Nucl. Phys. A 740 (2004) 
195. 

FP6 Euratom project EUROTRANS/NUDATRA, EC contract number 
FI6W-CT-2004-516520. 

A. Letourneau et al, Nucl. Phys. A 712 (2002) 133. 
http://www-nds.iaea.org/spallations 

F.E. Bertrand and R.W. Peelle, Phys. Rev. C 8 (1973) 1045. 

J. B. Natowitz, R. Wada, K. Hagel, T. Keutgen, M. Murray, A. Makeev, 

L. Qin, P. Smith, C. Hamilton, Phys. Rev. C 65 (2002) 034618. 

VI. Bogatin, VF Litvin, O.V Lozhkin, N.A. Perfilov, Yu.P. Yakovlev, 

Nucl. Phys. A 260 (1976) 446. 

B. N. Mekhedov, V. N. Mekhedov, Yademaya Fizika Vol.1 1, p.708 (1970) 
L. A. Currie, W. F. Libby, R. L. Wolfgang, Phys. Rev. 101 (1956) 1557. 
E. L. Fireman and J. Zhringer, Phys. Rev. 107, 1695 (1957) 

J. P. Alard, A. Baldit, R. Bran, J. P. Costilhes, J. Dhermain, J. Fargeix, L. 
Fraysse, J. Pellet, G. Roche, J. C. Tamain, II Nuovo Cimento, Vol. 30a, 
n2 (1975) 

K. Goebel, H. Schultes, J. Zhringer, Report CERN 64-12, 1964. 

S.G. Mashnik et al, J. Phys.: Conf. Ser. 41 (2006) 340. 

K. Ammon, I. Leya, B. Lavielle, E. Gilabert, J.-C. David, U. Helpers, R. 

Michel, Nucl. Instr. and Meth. B 266 (2008) 2. 

A. Guertin et al, Eur. Phys. J. A 23 (2005) 49. 

M. Enke et al., Nucl. Phys. A 657 (1999) 317. 

I. Leya, R. Wieler, J.-C. David, S. Leray, L. Donadille, J. Cugnon, R. 

Michel, Nucl. Instr. and Meth. B 229 (2005) 1. 

E. Le Gentil et al., Phys. Rev. Lett. 100 (2008) 022701. 

J.C. David and J.S. Hendricks, private communication. 



6 



